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TABLE 2 

Q{s), ^(7r) Fits of IRAS Data to Models 



Moments 



Fits to Full Sky 



^12(10) 



ail) 



'X.Tnin 



SCDM 
BCDM 
OCDM 

SCDM 
BCDM 
OCDM 

SCDM 
BCDM 
OCDM 

SCDM 
BCDM 
OCDM 
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Exponential 

Exponential 
Exponential 
Exponential 

Exponential 
Exponential 
Exponential 



1671?! 
1331^8 
167l|i 



109lf, 

-'■-'■-'■-40 

10711= 



Fits to Full Sky, Clusters Boosted 

^ L7A ^ LA 



1671^4 
15011? 



Fits to Northern Sky 

L Kfl 



109lt? 
111+35 

111_49 

96+" 



000 + 56 
■J-J-J- 183 



Fits to Southern Sky 

Q1 7+72 



o-, Y+72 
•^^ ' -85 
01 7 + 52 
^J- I _40 



SlStJjo 
180lf, 

160l=/, 



3i7ir, 

3001^3 
317111 



333 
317 



+44 

46 
+ 45 
38 



333lfo 



367lf, 



350 
333 



+47 

51 
+ 59 
37 



. + 50 
36 



317l_ 

3001II 

300111 



22.0 
21.3 
22.2 

15.8 
15.3 
16.0 

21.3 
20.5 
20.1 

13.1 
12.5 
12.3 



All errors given are 1 — it. 



26 Fisher et al. 



TABLE 1 

Fits of 100 Mock IRAS SCDM Catalogues to Models 



Moments 




{^12(4)> 


Error 


{^12(10)> 


Error 




Error 






SCDM 


Exponential 


323.6 


59.4 


175.6 


32.2 


441.1 


42.4 


32.3 


34.7 


SCDM 


Gaussian 


407 


78.5 


220.9 


42.6 


339.3 


22.2 


44.4 


57.7 


SCDM 


V = 1.5 


354.6 


66.6 


192.5 


36.1 


360.1 


28.6 


1 37.6 


43.4 


BCDM 


Exponential 


244.0 


55.2 


166.5 


37.7 


457.4 


50.0 


29.4 


37.3 


OCDM 


Exponential 


415.0 


71.6 


282.8 


48.8 


407.9 


40.8 


38.6 


41.3 


ACDM 


Exponential 


239.3 


57.5 


129.9 


31.2 


441.5 


44.1 


33.3 


38.9 


ACDM 


Spline 


248.3 


102.9 


134.7 


55.8 


393.9 


55.5 


56.7 


60.5 


Means are over 


the 100 SCDM Mock IRAS catalogues. 


^^12(4), V12 


(10), o-(l), and their 


errors (i.e., the standard deviations 


over the realizations) are listed 


in km s ^ . Th 


e two 


columns are means of the best 


solution. 


and means 


of the solution 


with the correct values of ■Ui2(4) 


and ij(l) imposed. The 


correct values of ■Ui2(4), 


^12(10), 


and ij(l) 


are 349.1, 189.5, and 420.3 



km s ^ , respectively. 
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the true cr(^r') is a strongly declining function of scale (much more so than shown in the OCDM simulation), and we have erred 
in modeling it as very shallow. It is difficult to see how this could be, particularly since the low density models of structure 
formation eschew bias in the galaxy distribution, and the measured cr(^r') in this model must exist at roughly the inferred value 
merely to satisfy the Cosmic Energy Equation and balance the potential energy stored in the galaxy fluctuations. 

Conversely, to accommodate a high density unbiased universe with our measured streaming, one must argue that cr(^r') 
grows more steeply with scale than has been modeled in our analysis. The strong covariance between cr(^r') and ■Ui2('') might 
then permit a larger estimate for the inferred streaming, which we presently estimate to be more than a factor of two smaller 
than expected in an Q = 1, unbiased universe (i.e. Ri 2a from /3 = 1). Our virialization procedure for cr(^r') on small scales 
insures it(4)/it(1) = 4} "'^ = 1.27. Given the strong constraint on ij(l), a larger asymptotic value of a[r > ro) would imply 
either a steeper gradient in cr(^r'), or a larger scale for the cutoff, or both. It seems dubious that sufficient power in the matter 
distribution on the scale of 5-10 /i~'^Mpc exists to generate this excess dispersion. Power on still larger scale (such as to provide 
for large amplitude flows on scales greater than 3000 km s"'') is ineffective at generating relative dispersion between pairs of 
galaxies on the scales under consideration here. 

7 SUMMARY 

The main conclusions of this paper can be summarized as follows: 

• Redshift distortions on small and large scales are observed unambiguously in the 1.2 Jy IRAS galaxy redshift survey. These 
provide a statistical measure of peculiar velocities, obtained purely from redshifts, with no independent measurements of 
distance. 

• The redshift correlation function ^{tp,t) can be modeled adequately on scales ^ 15 /i~''Mpc by a convolution of the real 
space correlations with the relative velocity distribution function. We can infer information about the first two moments 
of this velocity distribution. The velocity distribution function itself can be taken to be an isotropic exponential; skewness 
and anisotropy of the distribution function create distortions of opposite sign and roughly cancel one another. Much larger 
redshift surveys, such as the Sloan Digital Sky Survey (Gunn & Knapp 1993) will have the signal-to-noise ratio to allow 
us to look for these more subtle physical effects. This survey will cover a much greater volume than that of the IRAS 
survey, greatly decreasing the effect of individual structures in ^('•p,7r), and thus allowing a much cleaner measure of the 
redshift distortions. Modeling such higher quality data will require JV-body simulations which reproduce the correlation 
function of the survey more closely than the simulations used in this paper, eliminating the need to rescale the functional 
forms of the velocity moments, and reducing the possibility of systematic errors in the covariance matrix. 

• The mean relative streaming of IRAS galaxies on a scale of 4/i~'^Mpc is Vi2 (4) = km s"^ while the pairwise 
velocity dispersion at lh~^M.pc is ij(l) = 317^^g km . This streaming flow acts to cancel the Hubble flow; our 
measurements imply approximately half of the Hubble expansion of pairs is canceled on the scale of ro ~ 4 h~^M.Tpc. The 
mean streaming measurement has not previously been reported, while the pairwise velocity dispersion is consistent with 
previous estimates derived from optically selected redshift catalogues. The results are robust to double-counting cluster 
cores, but the amplitude of the mean streaming has high covariance with the derived pair dispersion. Moreover, it is quite 
sensitive to the model assumed for the dispersion on small scales. 

• The Cosmic Virial Theorem applied to the derived IRAS velocity dispersion leads to an estimate of Q/b ~ 0.38 ± 0.24. 
The relatively large uncertainty in this estimate reflects the statistical uncertainties of a nearly divergent integral over the 
three point correlation function. However, systematic errors are expected to be even more severe, and largely invalidate 
the Cosmic Virial Theorem as an estimator of cosmological parameters. 

• Linear theory and the derived estimate of the streaming give an estimate of /3 = Q°'^ /b = 0.45^Q ^g (statistical errors). 
The derived amplitude of the streaming is sensitive to the form assumed for the velocity dispersion on small scales. 
Substantial improvement of the error bars on this estimate will require much larger redshift surveys. 

• Our inferred /3 is applicable to an intermediate scale, ~ 10 — 15/i~'^Mpc, and we have obtained results intermediate 
between the lower values determined from estimates of cluster masses Q/fe = 0.1, r < l/i~''Mpc), and those obtained on 
much larger scale (/3 = 1.28, 15 < ?• < 40 /i"^Mpc). If all th ese results are correct, they argue for strong scale dependence 
of galaxy biasing. 
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Figure 14. Constraints on n and b from tlie derived velocity dispersion and large scale streaming. The results for different models 
for the velocity moments are shown with different line types, a). Estimates of Q/b from the Cosmic Virial Theorem and the derived 
normalization to <7(r) for the IRAS full-sky data. Because we use the Cosmic Virial Theorem to modify the AT-body <7(r) on small scales, 
the results are independent of scale from 1 to 3.76 Mpc. As discussed in the text, systematic errors make this determination worthless. 
6). Estimates of /3 from Equation 21 and the derived normalization to ui2(r) for the IRAS full-sky data. Because we use linear theory 
to modify the AT-body vi2{r) on large scales, the results are independent of scale above ~ 7/i~^Mpc. 

high density universe in which the bias in the galaxy distribution decreases strongly with scale. Scale dependent bias operating 
on scales of less than 1 /i^^'Mpc, where baryonic cooling and dynamical friction are effective, would not be too surprising 
(White & Rees 1978; Cen & Ostriker 1992; Couchman and Carlberg 1992). However, extending this scale dependence to larger 
scales is difficult to understand. Only more precise estimates of ■Ui2(7'), derived from much larger redshift surveys, can improve 
our error bars substantially. 

As discussed above, our estimate of ■Ui2(7') will have large systematic error if our assumed form of it(7') is grossly in error. 
The velocity dispersion differs greatly from its linear theory limit on all scales that we probe, thus assuming linear theory 
alone will cause systematic errors in any derivation of Vi2. This interdependency of the infall streaming and the nonlinear 
pair dispersion is inevitable, and in our assessment effectively invalidates Hamilton's (1993a) analysis which is based solely 
on the effects of linear theory. Hamilton extends his analysis to extremely large scales, but the expected infall streaming is 
very small on these scales, and as Figure 6 demonstrates, on all scales for which one can hope to measure the distortions in 
^(7'j,,7r), the role of the non-linear component of it(7') cannot be neglected. 

In order to reconcile an apparent large value of ■Ui2(4) with a low density universe, it would be necessary to argue that 
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6 DISCUSSION 

In § 3 we described how measurement of the redshift distortions of the ${rp, tt) maps provides two independent estimators of 
the cosmic density parameter Q. Both of these estimates are sensitive to bias relating the galaxy and mass distributions, but 
their dependency is different in the two cases. 

In order to apply the Cosmic Virial Theorem, we need a value for Q which expresses the ratio between the second and 
third moments of the IRAS galaxy density field. Bouchet et al. (1993), using a counts in cells analysis, derived a value for the 
"normalized" skewness, 53 = 1.5 ± 0.5, which is related to Q by Q = JS3J3/J2 where J2 and J3 are dimensionless integrals 
defined by Peebles (1980). For the IRAS value of 7 = 1.66 this yields Q = 0.48 ±0.16. This is roughly half the value obtained 
for optically selected galaxy samples (e.g., Groth & Peebles 1977), and is a reflection of the relative undercounting of cluster 
cores in IRAS selected galaxy samples (Strauss et al. 1992a). Indeed, Bouchet et al. found that the derived value of 53 (and 
hence Q) is extremely sensitive to the way in which clusters are treated in the analysis, and that giving them extra weight 
increased Q by a factor of two. 

The J(7) integral which appears in Equation 24 can be evaluated numerically for the IRAS value of 7 = 1.66 to yield 4.28. 
On small scales, ^ ■U121 thus the latter term can be dropped from Equation 24. Because we have used the Cosmic Virial 
Theorem to force the shape of cr(^r) at small distances, our estimate for Q/fe is independent of r from 1 to 3.76 h~^M.Tpc. Thus 
we find Q/b = 0.38. What is the error in this quantity? Statistical errors arise in the determination of it, 53, 7, and To, which 
lead to an estimated 1 — it error in Q/b of ± 0.24. Large as this statistical error is, it is overshadowed by systematic effects. 
The derived value of Q from the IRAS sample is exquisitely sensitive to the treatment of clusters; boosting them causes Q 
to double (Bouchet et al. 1993). Deviations of the two and three point correlation functions from their assumed power law 
forms will lead to further systematic error. + The Cosmic Virial Theorem is derived assuming statistical equilibrium of close 
pairs, and requires taking the continuum limit (Peebles 1980), neither of which can be justified observationally. Finally, as 
discussed earlier, it is only with perfect linear bias that the factor 'b' within the Cosmic Virial Theorem holds, although linear 
bias becomes an oxymoron on the highly non-linear scales on which the Cosmic Virial Theorem is valid (cf., Szalay 1988). We 
attempted to apply the Cosmic Virial Theorem to the JV-body models themselves with little success; the correlation functions 
differed too strongly from power laws to get meaningful results. Thus we conclude that we are unable to get a meaningful 
estimate of cosmological parameters from the Cosmic Virial Theorem, and in fact do not believe that the assumptions which 
go into it (perfect power law correlation functions, perfect linear biasing, isotropy of the velocity dispersion, the continuum 
limit, and statistical equilibrium of clustering) have been properly validated for any extant data sample. 

Although we cannot use our measurement of it(1) to obtain a useful estimate of Q, its value does put strong constraints 
on cosmological models. Standard CDM seems to require strong biasing to be consistent with the small value observed (Davis 
et al. 1985; Gelb & Bertschinger 1993), although it has been argued that if galaxies have systematically smaller velocities than 
do dark matter particles, the model can be saved (Couchman & Carlberg 1992). Recently, models with a mixture of hot and 
cold dark matter have become popular; these hybrid models better match the large scale constraints and at the same time 
have somewhat lower amplitude velocity fields on smaller scales, but it is not yet clear if their pair dispersion is an adequate 
match to the observations (Klypin et al. 1993; Davis, Summers, & Schlegel 1992). 

A more reliable estimate of cosmological parameters uses our measurement of ■Ui2(7') on large scales (Equation 21). Figure 4 
above showed that given the correct amplitude for ■Ui2(''), this method indeed yields the correct value of /3. Figure 14 shows 
the derived value of /3 as a function of scale from the normalization of V12 and our various models for the run of V12 with 
distance. Because we have scaled ■Ui2('') to agree with the linear theory prediction given the IRAS correlation function, the 
curves must asymptote to a constant value on large scales. However, the fact that the curves using different models largely 
agree is significant, and gives us confidence in the reality of the result. The curves give /3 = 0.45^0 ~ "')! the statistical 

error on this is dominated by that of ■Ui2(4). We do not include the possible systematic error due to uncertainty in the shape 
of it(7') on small scales (cf., the discussion at the end of the previous section). 

This measurement of /3 is on a length scale of 5-10/i"^Mpc, and is consistent either with the IRAS galaxies having a 
bias b = 2.2^q | in an Einstein-de Sitter Universe, or with the IRAS galaxies being an unbiased tracer in an Q = 0.26^q '^^ 
universe. From Table 2, we note that the Southern and Northern hemispheres alone yield estimates of ■Ui2(4) considerably 
higher than found for the full sky analysis. These estimates are consistent with /3 = 1, but with considerably larger error than 
for the full sky sample. Although all subsamples are consistent with each other within 1.5it, they span an enormous range of 
derived /3, indicating that the present redshift catalogs are too small for a definitive application of this statistic. 

Note that the allowed range of /3 lies between the typical values found on cluster and smaller scales from virial arguments, 
Q/b ~ 0.1 (or P ~ 0.25 if = 1), (Faber & Gallagher 1979; White 1990 and references therein) and the best fit value of 
P = 1.28 ± 0.3 (If) given by the POTENT-IRAS comparison on a scale of 15-40 /i"^Mpc (Dekel et al. 1993; see also Kaiser et 
al. 1991; Strauss et al. 1992c; Fisher, Scharf, & Lahav 1993). If all these observations are correct, then we apparently live in a 

+ The Cosmic Virial Theorem can be generalized for any form of the two and three-point correlation functions, but this would require 
knowledge not only of , but the full form of r, x — r), which has not been calculated for IRAS galaxies. 
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a) IRAS - Full Sample b) Best Fit SCDM Model 




Tp (h"^Mpc) 



Figure 13. a). The IRAS 1.2 Jy ti'^pt foi" the full sample as in panel a of Figure 1. b). The best fit SCDM model corresponding to the 
parameters in Equation 28 with an isotropic exponential velocity distribution function, c.) The difference between the IRAS and SCDM 
model. Contour levels are the same as in Figure 1. For clarity, the maps have been twice smoothed by a 1-2-1 boxcar in each direction. 
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Figure 12. Contours at larger than the minimum value (marked with a star) by 1.0, 2.3, and 4.0. These correspond to the 1 — a 
confidence level for the parameters taken separately, the 68% joint confidence interval, and the 2 — a confidence level for the parameters 
taken separately. These are based on the isotropic exponential fit to the whole sky IRAS data, using the SCDM model curves. The dots 
indicate the grid of values at which models were calculated. 
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An example of the constraint in parameter space (it(1), ■Ui2(4)) from the models is given in Figure 12, for fits to the 
isotropic exponential model using the SCDM velocity moments. The innermost contour describes the 1 — it error limits for 
each parameter taken separately (giving rise to the errors listed in Table 2). The second contour encloses 68% of the likelihood, 
while the third contour describes the 2 — it error limits for each parameter. There is strong covariance between the parameters, 
although it(1) is better constrained than is ■Ui2(4). 

Figure 1 shows that giving extra weight to clusters in the IRAS sample has only modest effects on ^('•p,7r). This is borne 
out by the model fits; the boosted sample shows slightly higher velocity dispersion than the unboosted case. This is reassuring; 
because the velocity dispersion is a pair-weighted statistic, we worried that our results would be quite sensitive to the way 
in which clusters are treated. It is possible, of course, that the galaxies in the very core of rich clusters, where the velocity 
dispersion is highest, are missing altogether from the IRAS sample, in which case no amount of cluster boosting would bring 
them back. However, Strauss et al. (1992a) found that the Finger of God of the Virgo cluster in IRAS and optically-selected 
galaxies had roughly the same extent in the redshift direction, arguing against this possibility. 

The velocity dispersion found in the South subsample is in good agreement with that for the full sky, while that for the 
North is somewhat larger, as we saw already from the ^('•p,7r) maps themselves in Figure 1. However, in both subsamples, the 
estimate of ■Ui2(4) is much higher than that in the full sky, even though the (5(s) statistics are in good agreement between the 
different subsamples (Figure 2); rather than choosing the compromise solution of the full sky that goes between the first six 
points of Figure 11a, the best-fit solution for the North and South subsamples goes through the second set of three points (like 
the lower dot-dashed curve of Figure 11a). However, the errors in ■Ui2(4) (at least in the North) are enormous, and include the 
value found for the full sky. From Figure 2 it is clear that there is appreciable noise in the estimates of (5(s) and ^(7r) even in 
the full sample; cutting the sample in half only amplifies the noise. Moreover, the elements of the covariance matrix for the 
half-sky mock catalogs are considerably larger simply because of counting statistics. These two effects combine to make the 
contours for the velocity moments in the Northern and Southern subsamples very broad, and the derived normalizations 
of ■Ui2(7') and (5(s) subject to large fluctuations. 

We conclude that the IRAS data are best fit by 

it(1) = 317I49 km , 

i;i2(4) = I67I7I km s"^ , (28) 
i;i2(10) = 109^47 km s"^ 

The velocity dispersion estimate is consistent with the value derived for the CfAl survey, 340 ± 40 km s"^ (D avis & Peebles 
1983), again implying that the undercounting of cluster cores by IRAS galaxies is not a serious issue. The streaming result 
has not previously been reported. 

The SCDM model corresponding to the parameters in Equation 28 is shown in Figure 136. This should be compared 
with the IRAS ^('•p, tt) from the full sample in Figure 13a (copied from Figure la). The difference of the two is shown as 
Figure 13c. The model is not perfect, particularly on small scales, but does provide a reasonable fit to the distortions on scales 
> 3 h-^Mpc. 

We have assumed the virialized form for it(7') on small scales for all analyses of the IRAS data. It is sobering to see the 
effect of dropping this assumption. We took two extreme models; in the first, we used the velocity moments of the SCDM 
model in which no smoothing on small scales had been done; so that it(7') rises on small scales. In this case, the oblateness 
of the ^('•p,7r) contours at large r is largely driven by velocity dispersion rather than infall (compare Figure 7), and the data 
are consistent with no infall at all (i.e., ■Ui2(4) = 0)! Alternatively, we use the SCDM model for ■Ui2(''), but for it(7') we use a 
form that grows more slowly than the virialized form: 

^ >~\a{l){10f'^ ifr>10fe-^Mpc. ^ ' 

In this case, the best fit values of ■Ui2(4) is now 417km s"'', more than twice the previous value. In both of these models, 
the best value of it(1) was essentially identical with that found above; moreover, all models gave comparable best values. 
We have no a priori reason to believe that it(7') differs substantially from a virialized form on small scales (indeed, Davis & 
Peebles (1983) present evidence that it(7') is a slowly rising function of ?• in small scales, close to the virialized form), but if it 
does, our estimate of ■Ui2(4) could be in error by a large factor. Our estimate of it(1) is quite robust to such uncertainties. 
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Figure 11. a). Plot of Q{s^ for the IRAS whole sky sample. The error bars are the standard deviation over the 100 SCDM mock IRAS 
catalogues, and are thus the same as are shown in Figure 9. The curves are the best fit normalizations given in Table 2 for various 
models. In addition, we show the model curves when the normalization of vi2 is moved up and down by its one-a errors (in the SCDM 
exponential fit). 6). Same as in panel a, for the ^(tt) statistic. 



modified to be in close agreement with one another (compare the left and right panels of Figure 3). Indeed, Table 2 shows 
nearly identical results for the normalizations derived from the various models. 

The two flanking dot-dashed curves in Figure 11 give the results of the SCDM model when the best-fit value of ■Ui2(4) is 
raised and lowered by its 1-a errors. The best-fit model cannot go through both the first set and second set of three points in 
the (5(s) plot, and strikes a compromise between the two. When ■Ui2(4) is raised and lowered by its errors, one obtains curves 
that go through one or the other of these sets. 

Table 2 shows that the minimum values of for our fits agree with each other quite well. This means that we cannot 
distinguish between the different models for the velocity moments, at least with the modifications made to them to bring 
them into closer agreement with the expected forms for the IRAS galaxies. However, the minimum values of are somewhat 
lower than we found with the mock IRAS catalogs, even when we used exactly the right form for the correlation function 
and velocity moments in the models. Indeed, Figure 10 shows that only 13% of the mock IRAS catalogues give a x^ ^ 22.0. 
The SCDM model has a correlation function which is substantially steeper than that of IRAS galaxies, which means that the 
covariance matrix may have overestimated the errors of the observed correlations on small scales. Although this is unlikely to 
strongly bias the best-fit models (it is hard to see how the model fits in Figure 11 could be improved) it will bias the derived 
value of x^ down, causing us to overestimate the errors on derived parameters slightly. 
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Figure 10. Cumulative distribution of the best (solid histogram) and the with normalization fixed at the correct value (dotted 
histogram) from fits (using an exponential velocity distribution function and the SCDM velocity moments) to Q{s^ and ^(tt) of the 
100 mock IRAS catalogs. The light and heavy smooth curves are the expected x^ distributions with 27 and 32.3 degrees of freedom, 
respectively. 

also shows results of fitting the 100 mock catalogues to both a Gaussian (j/ = 2 in Equation 12) and hybrid = 1.5) velocity 
distribution function and to the velocity moments from the open (OCDM) and biased (BCDM) simulations. Also shown in 
Table 1 is the results for the SCDM model with both the skewness of the velocity distribution (/(wal'') given by the spline fit 
in Figure 5) and the anisotropy of the velocity dispersion taken into account; we refer to this model as ACDM. 

The mean values show that the exponential model is greatly favored over the other models for the velocity distribution 
function. This is because ^(7r) gives us such a clear signature of what the correct functional form of /(raa) must be. Indeed, the 
exponential model is preferred over the spline fit of Figure 5, presumably for the reasons discussed in reference to that figure. 
The mean of the derived velocity dispersion at 1 /i~''Mpc is quite insensitive to the model used, giving results close to the right 
answer. However, the derived streaming at 4/i~'"Mpc does depend sensitively on both the shape of the velocity distribution 
function, and more ominously, on the form of the velocity moments used. Figure 9 shows that although the normalizations 
of the velocity distribution functions are quite different in the different cases, the model fits to (5(s) and ^(7r) are almost 
indistinguishable from one another, and are in good agreement with the data. In addition, when we fit models to the IRAS 
data, we will apply the modifications to the velocity moments discussed in § 4.1 above, bringing them into closer agreement 
with each other and thus reducing the discrepancy in normalization between them. Finally, the linear theory expression for 
Vi2 of course only holds on large scales; note that the normalization found for the BCDM and SCDM models (the two with 
the best x^) agree much better at 10 /i~''Mpc than at 4/i~''Mpc. Figure 4 shows that with the correct normalization, Vi2 does 
indeed yield a correct estimate of /3, giving us confidence that we will derive an unbiased estimate of /3 from the IRAS data 
themselves. 



5 MODELING THE /i?^5 ^(7-p,7r) 

Having developed statistical procedures that give meaningful answers when tested on JV-body models, we are ready to apply 
them to the IRAS data themselves. We have fit series of models to the whole sky IRAS data, the IRAS data with clusters 
boosted and to the Northern and Southern subsamples, with results listed in Table 2. For each fit, the 1-a errors on the two 
normalizations are given, as calculated from the projection of the Ax^ = 1 contour. In addition, we give the value of ■Ui2(7') at 
?• = 10 h~^M.pc, a scale that is fully in the linear regime, as Figure 4 shows. Finally, we list Xmin: the value of x^ for the best 
fit model. The fits to the North and South subsamples used covariance matrices calculated from 100 mock IRAS catalogues 
covering half the sky. 

Figure 11 shows the ^(7r) and (5(s) curves for the full-sky sample, together with the expected 1-a error bar per point from 
the mock IRAS catalogues. The three almost indistinguishable model curves shown are from the best-fit normalization of the 
SCDM, BCDM, and OCDM models with an exponential velocity distribution function. The models agree with one another 
better than was found for the fits to the JV-body models, because the current fits use velocity moment shapes that have been 
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Figure 9. a). The solid points plot Q{s^ for the average of the 100 mock catalogues, while the error bars show the dispersion around this 
mean. The curves show the mean of the best fit models to each of these catalogues as given in Table 1, for various forms of the velocity 
distribution functions. 6). The same as in panel a, now for the ^(tt) statistic. 



we know the model for the velocity distribution function is not correct in detail, since we have ignored its skewness and 
the anisotropy of the velocity dispersion tensor. However, we saw above that the effects of these two refinements on ^('•p,7r) 
approximately cancel each other, and ignoring them does not introduce any strong bias in the derived parameters. In fact, 
our effort modeling them greatly worsens the fit. As mentioned above, this is probably due to our overestimation of the true 
skewness; we have not attempted to refine the these anisotropic models, and we drop them from further attention. Because 
of the relatively poor signal to noise ratio of an individual realization, the simple isotropic exponential model for the velocity 
distribution is more than adequate, and will serve as our preferred model for the available IRAS data. As much larger redshift 
surveys become available, it may become desirable to refine our model by including additional effects within the velocity 
distribution function. 

Thus far, we have shown that when we use the correct functional form of cr(^r') and ■Ui2('') for the simulation, we can 
recover their correct normalizations. Of course we do not know the functional form of the velocity moments for the IRAS 
data. Therefore our ability to model real data rests on the insensitivity of our model ^('•p,7r) to the exact functional form of 
cr^r) and ■Ui2(''). We saw in Figure 3 that the scale dependence of the velocity moments for different cosmological models was 
similar although not identical. As a test we have tried to recover the normalization of the velocity moments of the SCDM 
mock IRAS catalogs with various other choices for the model of the velocity moments and distribution functions. Table 1 
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where the i and j range over the individual points in the {jp, tt) plane at which i{Tp, tt) has been measured. Our first approach 
was to minimize the quantity with respect to the normalizations of the CT(r^ and ■Ui2('') curves, for a given shape of these 
curves and a given form of the velocity distribution function. There are, however, several problems with this direct approach. 
The covariance matrix is difficult to measure accurately, even with 100 realizations, with the consequence that the estimated 
M can be close to singular. Furthermore, the error distribution of the 100 individual realizations shows significant skewness 
over all of the (?•;;, tt) plane, while the statistic above is meaningful only when the error distribution is a multi-variate 
Gaussian. In addition, direct fits to the ^('•p, tt) map are necessarily very sensitive to the amplitude and slope of ^(r), and to 
small deviations from a power law, while we are most interested in fitting the distortions of the ^('•p,7r) contours. Presumably 
because of these problems, we are unable to recover the correct normalization to the velocity moments for the mock IRAS 
catalogs (in particular, our estimate of the normalization of ■U12 is biased low by ~ 20%). 

With this in mind, we carry out our fit not to the full ^('•p,7r) contour map, but to the two projections we defined above 
in §2.2. As we described there, ^(7r) is a statistic that is almost entirely determined by the amplitude of cr(^r'), while (5(s), 
although sensitive to cr(^r'), contains most of the available information on ■Ui2(r) (recall Figure 7). Moreover, both statistics 
are normalized in such a way as to minimize their sensitivity to ^(r) itself. We find the statistics of (5(s) and ^(7r) to be well 
behaved over the 100 mock IRAS catalogues, with skewness of the error distribution consistent with zero on all scales, giving 
us hope that fitting the analogue of Equation 27 to the data will give unbiased results. 

In fitting models to data (either the actual or mock IRAS catalogues), we evaluate ^(7r) every 1 /i"^Mpc from 0.5 to 
14.5 /i~'^Mpc, and evaluate Q{s) every l/i"^Mpc from 1.5 to 15.5/i"^Mpc. The offset arises because it is not possible to 
measure (5(s) with our binned data on smaller scale, and for convenience we want to keep the vector lengths identical. The 
covariance matrix of these 30 quantities is determined from the same 100 mock IRAS catalogues discussed above. The quantity 
X^ is determined using Equation 27 above with ^('•p,7r) replaced with the vector containing Q{s) and ^(7r). Actually, because 
of the normalization applied to ^(7r) (Equation 7), the covariance matrix is close to being singular and the calculation of the 
inverse of the covariance matrix in Equation 27 is unstable. We circumvent this difficulty by using using Principal Component 
Analysis, as described in the Appendix of Paper 1. We find linear combinations of the data points which have a diagonal 
covariance matrix. These linear combinations are simply the original data vector multiplied by the matrix whose columns 
are the eigenvectors of the original covariance matrix. Furthermore, the diagonal elements of the covariance matrix of the 
new data points are merely the eigenvalues of the original covariance matrix. The near singularity in the covariance matrix is 
manifested as one eigenvalue which is much smaller in magnitude than the others. We discard the linear combination of data 
points which gives rise to this eigenvalue, thereby making the fitting procedure robust. Rejection of the smallest eigenvalue 
reduces the number of degrees of freedom in the fit by one for a total of 27 (30 data points - 2 parameters - 1 rejected). 

As a prelude to fitting models to the actual IRAS data, we can check the validity of our procedure using the mock IRAS 
catalogues. We calculate a grid of x^ values as a function of the normalization of cr(^r) and Vi2{r), for each of the 100 mock 
IRAS catalogues. Here we use an exponential model for /(wair), which we saw above accurately reproduces the mean ^('•p,7r) 
of the mock catalogues, we use the ^(r) of the SCDM model (as measured from the JV-body model itself), and we use the 
shapes of cr(^r') and ■Ui2(r) of SCDM (the solid curves in Figure 3 above, without correction for the virial theorem on small 
scales). The results are shown in the first entry of Table 1. We express the normalizations of the velocity moment curves in 
terms of the velocity dispersion at 1 /i~''Mpc, ij(l), and the streaming at 4/i~''Mpc, ■Ui2(4). For each of the 100 realizations 
we minimized x^ to find the best amplitude for the velocity moments; the mean of the best fit parameters are 323.6 km 
and 441.1 km s"\ for a{l) and ■Ui2(4) respectively, within 7% of the correct answer (which we know a priori from the JV-body 
mode 1) of 349.1, 420.3 km s"\ The rms scatter from realization to realization of these normalizations are 59 and 42 km s , 
respectively, and the covariance between these quantities is small. The mean ^(7r) and (5(s) over the 100 mock catalogues are 
shown as the points in Figure 9, with error bars given by the standard deviation over these catalogues. The solid curve is the 
best fit isotropic exponential model. The other curves will be discussed below. 

Figure 10 shows the cumulative distribution of the derived values of x^ for the hundred mock IRAS catalogues, for the 
best fit for each model (solid histogram, mean value 32.3) and when cr(^r') and ■Ui2(r) are held at their correct values (dotted 
histogram, mean value 34.7). If our model is a correct description of the data, the solid histogram should agree with the x^ 
distribution for 27 degrees of freedom, which is shown as the light smooth curve. They are clearly in poor agreement; the true 
distribution is better fit (although not perfectly!) for 32.3 degrees of freedom (heavy smooth curve), implying that our model 
is not a perfect match to the data. We have repeated this experiment limiting the fits to the range 1 - 10/i"^Mpc in Q{s) 
and ^(7r), as the noise gets very large beyond this range; the excess in the mean x^ over the expected value remains, although 
the derived ■Ui2(4) and ij(l) are in agreement with those listed in Table 1. 

Thus our model is not a perfect descriptor of the JV-body data; our average x^ is larger than the number of degrees 
of freedom by five, meaning that there are physical effects which we have not taken into account. Nevertheless, we have 
demonstrated that minimizing x^ gives unbiased results for the cr(^r) and Vi2{r) normalization. This is reassuring, because 



Clustering in 1.2 Jy Survey 15 



Average of 100 Mock IRAS Catalogues Exponential - <SCDM> 

20 L I J I r I — I — I — I — I \ \ rr^ \ 1 I I rri i 20 nn \ ^ \ \ \ rrn i i i i 




15 20 



o 




Fp (h"' Mpc) 



Fp (h"' Mpc) 



20 



SCDM, Skewness and Anisotropy 

I I I I I I I I I ^ i I i 

c) 



SCDM, Skew - Exponential 




5 10 15 

F (h"' Mpc) 



20 




5 10 15 

F {h-^ Mpc) 



20 



Figure 8. a). Average ^(t'pjTt) map for the 100 mock IRAS SCDM catalogues. Compare with the model in Figure 6a. 6). The difference 
between the model in 6a and the average of the realizations in panel a). All substantial differences are close to the origin, c). Model 
^(rp,7r) obtained by using a spline fit to the unbiased SCDM velocity distribution function, /(u^alr) (the heavy solid curve in Figure 56). 
In addition, the anisotropy of the velocity dispersion tensor was explicitly taken into account. Compare with Figures 6a and panel a. cf). 
The difference between the anisotropic skew model in panel b and the isotropic exponential model in Figure 6a. 



This tends to increase the prolateness of ^(rp.Tr) on small and intermediate scales. Thus these two effects work in opposite 
senses, and when both are included in the modeling, we get the model ${rp, tt) shown in Figure 8c. This model is quite close to 
the simpler model in Figure 6a. The difference between these two models (Figure 8d) is somewhat greater than that between 
the isotropic exponential model and the mean of the 100 mock IRAS catalogues; that is, the anisotropic skew model is a worse 
fit, even though it includes more relevant physical effects. We suspect that this is because this model overestimates the true 
skewness for the reasons discussed in reference to Figure 5. Although we do not show it here, adding either the refinement of 
skewness or anisotropy alone to the model causes a much larger change in ^('•p,7r) than putting both in; that is, they largely 
cancel one another's affects. 

As we show next, the exponential model with isotropic velocity dispersion does give unbiased results for the amplitudes 
of the velocity moments, and thus we will use it as our primary model to compare with the data. 

4.3 Fitting Models to the Mock Catalogues 

Given the covariance matrix M, we define the of the difference between the data and a given model: 
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Figure 7. a). The statistic Q{s^ derived from the isotropic exponential model. The solid points are with the correct normalization for 
a'(r), while the open circles and stars are for <7(r) halved and doubled, respectively. These correspond to the three models of Figure 6. 6). 
The statistic ^(tt) derived from the three models in Figure 6 and panel a. c). The statistic Q{s^ derived from the isotropic exponential 
model. The solid points are with the correct normalization for ui2(r), while the open circles and stars are for ui2(r) halved and doubled, 
respectively, cf). The statistic ^(tt) derived from the three models in panel c. 



ensemble of realizations, as described in the Appendix to Paper 1. Of course, the volume of the JV-body simulation itself is 
finite, and the same structure will appear in many of the realizations. However, it will appear with a different orientation in 
each one, and thus its coupling with the redshift space distortions will differ as well. Thus we feel justified in claiming the 
individual realizations to be independent. 

Before using these realizations to determine the covariance matrix, we can compare them with models. Figure 6a above 
showed a model for ^('•p,7r), using velocity moments correctly normalized for the SCDM model. Figure 8a shows the average 
over the 100 realizations of ^(t-j,, tt). Figures 6a and 8a should agree if the model assumptions that went into Figure 6a were 
correct; indeed. Figure 86 shows that the differences between these two are substantial only in the region immediately around 
the origin, where the correlations are strong. However, there are two physical effects we have ignored in our model. First, the 
true velocity distribution function differs from the exponential assumed above, as Figure 5 shows us; there is non-negligible 
skewness, whose effect will be to increase the oblateness in ^('•p,7r) on intermediate scales (0.5 ro ^ 2ro). The approximate 
universal velocity distribution function given in Figure 56 can be included directly in the model. Second, the velocity dispersion 
tensor is anisotropic; that is, the second central moment of the longitudinal component of velocity is larger than 2~^^^ times 
that for the transverse component. For the JV-body models, we can measure this anisotropy, and put it directly into the model. 
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Figure 6. Model ^(t'p, tt) derived using the exponential velocity distribution and the velocity moments of the unbiased SCDM simulation. 
Panel a) shows the resulting model using the velocity moments has they appear in Figure 3. In panel 6) the amplitude of the dispersion 
has been doubled while in panel c) the amplitude has been reduced by a factor of two. 



4.2 The Errors in ^(r"p,7r) and Mock IRAS Catalogues 

We now have all the ingredients necessary to make models to compare with the observed ^('•p,7r): the real space correlation 
function (Equation 4), a family of models for the velocity distribution (Equation 11), and an empirical determination from 
JV-body simulations of the scaling of the first and second moments of the velocity distribution with separation (Figure 3). 
In order to fit these models to the data, however, we need to assess the statistical errors in ^('•p,7r). In Paper 1, we used 
bootstrap methods to assess the covariance matrix of the redshift space correlation function, although it became clear there 
that the bootstrap technique systematically overestimates the correlation function errors (c.f.. Mo, Jing, & Borner 1992; 
Press et al. 1992). In addition. Figure 1 dramatically shows that the number of coherent structures included in the IRAS 
volume is small enough to leave their signatures in ^('•p,7r). These same structures will be present in all bootstrap resampling 
realizations, and thus the derived error will not reflect the uncertainty introduced by these structures. With this in mind, we 
take a different approach to defining the covariance matrix of ^('•p,7r). From the SCDM JV-body realization, we create 100 
realizations of the IRAS redshift survey, in the manner described in Paper 1. In brief, for each one, we choose an JV-body 
point with peculiar velocity, local density, and local shear similar to those of the Local Group, and select objects around it 
with number density as a function of distance as given by the IRAS selection function. The ^('•p,7r) is calculated for each 
realization, and the covariance matrix M of the determination of ^ at a grid of points in the ('•p,7r) plane is found from the 
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Figure 5. The distribution function, \t), of the longitudinal component of the relative velocity measured from pairs of points in the 

SCDM simulation, a). \t) as a function of the longitudinal component of the velocity difference, . The various curves correspond 

to /(u;|||r) at separations r corresponding to the clustering amplitude ^(r) = 50, 10, 1, 0.5, and 0.1 as labeled in the figure itself. 6). 
The same curves shown in a) scaled by their mean and variance a. The heavy solid curve denotes the mean of the distributions at 

different separations. Note the linear scale on the abscissa. 



velocity moments as plotted on the left side of Figure 3. In Figure 66, we have doubled cr(^r') from that in Figure 3a, and in 
Figure 6c we have halved cr(^r') from that in Figure 3a. Figure 6a shows considerable oblateness on large scale and moderate 
prolateness on small scale. Figure 66 is quite prolate on all scales, while Figure 6c shows extreme oblateness, especially on 
large scale. Figure 7a shows the derived for these three maps. The for the low dispersion model goes very negative 

(oblate) on large scales, while the high dispersion model has positive Q (prolate) on nearly all scales. Thus although in linear 
theory (5(s) is primarily a measure of the first moment Vi2 of the velocity distribution function, it is in fact strongly dependent 
on the amplitude of the velocity dispersion, which on small scales, (?• /S, 10 — 15 /i~''Mpc), is dominated by nonlinear effects. 
Therefore estimates of /3 from redshift distortions in ${rp, tt), which have a measurable signal only for modest scale separations, 
must necessarily model the dispersion accurately on scales where linear theory is not applicable. 

Figure 76 shows the ^(7r) plots for these three models; the principal effect of changing the normalization of a is to put 
an offset into (5(s), and to change the slope of ^(7r). Figures 7c and d show (5(s) and ^(7r) as the normalization of ■Ui2('') is 
halved and doubled. It affects the slope of (5(s), but ^(7r) is almost independent of the amplitude of ■Ui2(7'), which means that 
it gives a clean measure of the velocity dispersion, as it was designed to do. 
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Figure 4. Plots of j;i2(t-)(1 + ^(T-))/((2/3)fl"c,T- ^(t-)) versus r for the SCDM, BCDM, and OCDM models. In linear theory (large scales) 
this curve should take on the constant value /3 = /(n)/6. 



We consider the SCDM model in further detail to address this question. Figure 5a shows the longitudinal component of the 
relative velocity distribution function of pairs /(ranlr) extracted from this JV-body model, for various separations r, labeled 
by the amplitude of the clustering strength, ^(t'). This is part of the underlying distribution from which Vi2{r') and cr(^r') 
of Figure 3 were computed. The distribution is quite skew-negative on all scales. The distribution function appears to have 
exponential tails, but has a softer than exponential core. If these curves are functions only of their first and second moments, 
then they will all coincide when scaled to a common mean and variance (Equation 11). Figure 56 applies this scaling, and 
indeed shows that they are close to being self-similar. The heavy solid line is a spline fit to the mean of these curves; this mean, 
together with the measured ■Ui2('') and cr(^r'), gives a reasonable approximation to the true /(ranlr) for the JV-body models. The 
transverse component of this distribution, f{w±\r'), cannot, by symmetry, be skew, and the projected distribution, flwslr), 
is a combination of the two. Thus models using this spline form for /(wa) overestimate the skewness; our modeling below 
ignores this subtlety. 

Our aim is to measure /3 from the measured amplitude of ■Ui2(7') using the linear theory prediction of Equation 21 above. 
Unfortunately, on all scales probed by the present data set, the effect of ■Ui2(7') on ^('•p, tt) is strongly counteracted by that of 
cr^r). That is, we cannot determine ■Ui2('') without also determining cr(^r'), so we must simultaneously model both the linear and 
non-linear aspects of the problem. This is illustrated in Figures 6a, b, and c, which show models of ^(7'p,7r) from Equation 12 
using an exponential velocity distribution and the SCDM model for ^(r) and the velocity moments. Figure 6a uses the SCDM 
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Figure 3. Velocity moments for tlie AT-body simulations. The different line types give results for the different simulations used, as is 
indicated in the figure itself, a). The one dimensional relative velocity dispersion, 6) The mean pair-wise velocity, ui2(r). The 

diagonal straight line in the Hubble line. c). The quantity after the modifications to match the IRAS data on small and large 

scales, as described in the text. cf). The quantity ui2(r), after the modifications to match the IRAS data. 

galaxies forming pairs on smaller scales are taken from high density, high dispersion, environments. This is why a peaks at 
a separation characteristic of clusters, ~ 2 — 3/i~'^Mpc. On larger scales the dispersion declines to its asymptotic value of 
(2/3)'"''^ of the rms value of the 3-d velocity dispersion. 

The behavior of ■Ui2('') is characteristic of that seen in all JV-body experiments (e.g. Efstathiou et al. 1988): on large 
scales the mean infall streaming does not fully compensate for the Hubble expansion (the straight line of Figure 3b and d) but 
does retard the physical expansion of the pairs, while on smaller scales, r ^ro, the mean infall streaming overcompensates 
the Hubble flow, indicating that, on average, pairs on this scale are physically moving together. 

The similarity of the shapes of these curves to one another means that we will not be able to distinguish observationally 
between them; we will only be sensitive to the amplitude. In particular, we will not readily be able to distinguish biased 
and unbiased models with the same values of /3 = Q°'^/fe from the velocity moments themselves. This is illustrated directly 
in Figure 4, which shows the quantity 3/2'Ui2(7')(l + ^Ir)) / [H o r ^[r)) versus r, which, according to Equation 21, should 
approach the constant value of /3 on large scales. Indeed, in each case, the curves approach the correct value of /3 on scales 
above 7/i"^Mpc or so. This is small enough that we have a reasonable chance of detecting the linear theory limit in Vi2 in 
the IRAS data itself. 

Thus far we have considered the moments of the velocity distribution function. What is the best model for its shape? 
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4.1 Velocity Moments 

We have derived velocity moments from three variants on the Cold Dark Matter (CDM) model of structure formation. We 
will see that the shape (but not the normalization) of the moments as a function of r varies little from model to model, which 
suggests that they are generally applicable (Efstathiou et al. 1988) and can serve as templates for our modeling. 

The first simulation employed is that of a standard CDM model (Q = 1, Qh = 0.5 and Qa = 0; hereafter SCDM); 
the details of this simulation are given in Frenk et al. (1990). We compute the correlation function and velocity moments of 
both unbiased and biased subsets of the particles. The unbiased statistics are derived from a random fraction of the 262,144 
simulation particles at the time when the correlation amplitude of the points is close to that of the IRAS galaxies; we designate 
this epoch as redshift z = 0. The correlation statistics in the biased model (hereafter BCDM) are computed by weighting each 
simulation particle by the expected number of peaks per particle, using the "peak-background split" procedure of Bardeen et 
al. (1986). Galaxies are associated with 1.5 it peaks at an output time corresponding to z = 0.6, using the methods described 
in White et al. (1987). The correlation functions ^(r) of the BCDM and SCDM models are very similar, so the effective bias 
of the BCDM model is = (1 + z) = 1.6. We also examined the velocity moments of a simulation of an unbiased open CDM 
model [Qh = 0.14, Qa = 0; hereafter OCDM), kindly provided to us by G. Kauffmann from the model "M" simulation of 
Kauffmann & White (1992). We will measure the moments of the velocity distribution function from all three JV-body models, 
but we measure the shape of F and the covariance matrix of (5(s) and ^(7r), and test our fitting method, using the SCDM 
model alone. 

A common problem with coUisionless JV-body simulations is that their small-scale velocity dispersions are much higher 
than indicated by observations, including those presented in this paper (e.g., Davis et al. 1985; Gelb & Bertschinger 1993). 
This is often taken as an argument against unbiased Q = 1 models, although including the non-gravitational physics of 
baryonic material might substantially decrease the velocity dispersions of pairs (however, see Cen & Ostriker 1992). In any 
case, we wish to produce models that closely mimic the data, and thus we convolve the velocity field of the points in the 
simulations with a Gaussian window with a width of 0.25 /i"^Mpc. That is, our aim is not to test the underlying assumptions 
that go into generating the JV-body models, but simply to obtain generic fitting functions for the velocity moments. 

In applying these moments to the model fitting of IRAS data, there are three modifications we would like to make, in 
order that they more closely match what we expect given the clustering properties of IRAS galaxies, and to bring the curves 
into closer agreement with one another. First, the correlation functions of the IRAS galaxies and the JV-body particles are not 
the same; in particular, they have slightly different correlation lengths (7-0=3.35, 4.3, and 5.0/i"^Mpc for the SCDM, BCDM, 
and OCDM models respectively, while the IRAS correlation length is To =3.76 h~^M.pc). Thus we scale both the r coordinate 
and amplitude of the velocity moments by the ratio of the correlation lengths. In addition, the correlation functions of the 
JV-body models are somewhat steeper than that of IRAS galaxies, which will influence the small scale behavior of cr(^r'). The 
Cosmic Virial Theorem (Equation 24) tells us that the dispersion of bound objects scales like 

o-(x(5M/7-)'''= oc (7-=^(7-))''" ocT-'-'^''^ oc/ " , (25) 

where the last proportionality holds for the IRAS value of 7 = 1.66 (Paper I). We therefore replace the small scale velocity 
dispersion calculated from the JV-body models on scales r <C To = 3.76/i~'^Mpc with the scaling given in Equation 25; we 
furthermore hold a constant for r < l/i~''Mpc, as we find we get more stable estimates of the normalization of ■Ui2(7') when 
we do this. 

Finally, on large scales. Equation 21 shows us the scaling of V12 expected in linear theory. Because the JV-body models 
and the IRAS galaxies have different correlation functions, the slope of ■Ui2(7') will differ. We thus modify the JV-body ■Ui2('') 
as follows: 



^^12('') Vi2{t) 



1 + {^Y [2Pri{r)] I ^2(0(3 - 7)(1 + i{r))] 



(26) 



where ^ and 7 are appropriate for the IRAS galaxies (Equation 4), /3 is the value appropriate for the JV-body model in 
question, Tc = 5/i~''Mpc, and the form (^) is an arbitrary but convenient form to force a sharp transition to the correct 
linear theory slope on scales above Tc. 

Figure 3 plots ■Ui2('') and ct{t) as measured directly from the three simulations. The Gaussian smoothing of the velocity 
field is evident in the suppression of it on small scales ^ 1 /i~'"Mpc. The left two panels show the velocity moments as measured 
from the simulations directly, while the right panels have had the IRAS rescaling of lengths and the further modifications of 
Equations 25 and 26 applied. With these modifications, although the three curves in each case have very different amplitudes, 
the shapes are similar. 

The unbiased SCDM model has the largest amplitude velocity field, as expected. Note that it(7') is essentially flat on 
large scales for all models, and indeed, when the IRAS virialization condition is applied on small scales, it(7') remains constant 
to first order on all scales. Note that it(7') is the rms velocity dispersion about the mean for all pairs of a given separation. 
Most galaxies forming pairs of large separation are in low and medium density environments, while a higher proportion of 
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What is the signature of large-scale streaming in ^('•p,7r)? If the dispersion is a slowly varying function of scale (a good 
approximation on large scales, as we shall see in § 4.1) and ^ <C 1, then Equation 12 can be expanded in vi^ij)! HqT ^° S^^^ 



Hos 



, , 2d \vi2{s) 

as s 



+ 



V12 
HoS 



(22) 



where = r^ + T^, fj, = t/s, and ^(s) is the direction-averaged redshift space correlation function (cf., Peebles 1980, equation 
76.9). In the linear regime where ■Ui2('') is given by Equation 20, Equation 22 can be written as 



(23) 



Thus the dominant effect of the streaming is to introduce a quadrupole moment to ^('•p,7r) (which is the motivation for 
defining (5(s) as we did in the previous section). Although our expansion of ^('•p,7r) to first order in vi2{t)I HqT is not 
equivalent to a first order expansion in the density field 6, Equation 23 agrees to first order in /3 with the result of Hamilton 
(1992) (his linear theory expansion also includes a hexadecapole term). In particular, the angle average of Equation 23 gives 
{ajpt"^)) IJ- = (1 + (2/3)/3) ^('•), which is Kaiser's (1987) result to first order in /3 (Equation 1). The actual linear theory limit 
of Equation 12 differs from Equation 22 because we have ignored the scale dependence of the linear theory dispersion here. 
However, our modeling of ^(7'j,,7r) below takes this scale dependence into account explicitly, meaning that our results are 
consistent with linear theory on large scales (cf.. Figure 4). 

Equation 21 points us towards our first goal; if our model fits to ^('•p,7r) can give us ■Ui2(''), we can put constraints on 
p. Another constraint on Q and h comes from the second moment, ct{t). In the continuum limit, the statistical equilibrium of 
clustering on small scales implies that the tidal gravitational acceleration acting on pairs of galaxies is balanced by the kinetic 
pressure of their relative dispersion (Peebles 1980). Modeling the two-point correlation function as a power law of index 7 
and the three-point correlation function as a power law of index 27 yields the Cosmic Virial Theorem: 



S^(7-) - ■Ui2(7-) 



9QQ(g„7-)^^(7-)J(7) 
4fe(7-l)(2-7)(4-7) 



(24) 



where Q (not to be confused with (5(s) defined in the previous section!) expresses the proportionality between the two- 
and three-point correlation functions and J(7) is the dimensionless integral of Peebles (1980, Equation 75.13). The poles in 
Equation 24 for 7 = 1 and 7 = 2 attest to the sensitivity of the result to the power law form of the three point correlations. 
This expression is derived assuming that the velocity dispersion tensor is isotropic, and thus the three-dimensional velocity 
dispersion S is Z^^'^a. The dependence of the Cosmic Virial Theorem on biasing is ambiguous. As discussed by Bartlett & 
Blanchard (1993), the velocity dispersion of galaxies depends on the ratio of the galaxy-galaxy-mass three point correlation 
function to the galaxy-galaxy two-point correlation function. If we assume linear bias. Equation 24 follows. However, the 
Cosmic Virial Theorem is valid on small scales where the density contrasts are much larger than unity and linear bias is 
unlikely to be applicable. Furthermore, the bias appearing in Equation 21 is on large scales where ^ <C 1, and there is 
no guarantee that it has any direct relation to the small-scale bias. Indeed, Bartlett & Blanchard (1993) show that the 
interpretation of the quantity Q/fe which we solve for using Equation 24 in fact depends sensitively on the biasing model 
assumed. However, we show in § 6 that the combination of statistical and systematic uncertainties entering into the Cosmic 
Virial Theorem make it almost useless as an estimator of cosmological parameters, and thus the ambiguity of the role of 
biasing becomes a moot issue. 



4 CALIBRATION USING JV-BODY SIMULATIONS 

In order to model the observed ^('•p,7r) (or equivalently, its projections defined in § 2.2) from Equation 12, several ingredients 
are needed. The real space correlation function oi IRAS galaxies is known (Equation 4), but the velocity distribution function 
as a function of separation is not. Our approach is to measure the distribution function and its moments from a series of 
JV-body models, and use this directly in our modeling, varying only the form of the distribution function and the amplitude 
of its moments. This approach is valid to the extent to which the derived best amplitudes are insensitive to the details of 
the velocity moments adopted. The JV-body simulations themselves, and their velocity moments, are presented in § 4.1. We 
also extract mock IRAS redshift surveys from one of the JV-body simulations, from which we measure (5(s) and ^(7r), and 
thus quantify the covariance matrix between the observational quantities (§ 4.2). This allows us to fit models to the data by 
minimizing x^; in § 4.3, we test this by applying it to the mock redshift surveys themselves, and demonstrate that we recover 
unbiased estimates of the amplitudes of the velocity moments. We emphasize that we are not using these JV-body simulations 
to test the models that underlie the simulations; rather, we use them as a bridge over our ignorance of the true form of the 
velocity distribution function, for our modeling of ^(7'j,,7r) and its errors, and in order to test our fitting methods. 
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1980) and can be readily measured in JV-body simulations. Their qualitative behavior is insensitive to the power spectrum, 
as we will see in § 4.1 below. 

We parameterize the functional dependence of /(tosIt") on its first and second moments by writing 



f{W3\r) = —r^g 

cr[r) \ '^vi 



(11) 



In particular, consider the family of models in which g is an exponential of a power of its argument. Let the vector r separating 
two galaxies be decomposed into a transverse component Tp and a line-of-sight component y. Then = -k — y and Equation 9 
gives 



1 + ^('-p,t) = C 



f dy 
J <^ir) 



(1 + ^(7-))exp 



-V 



- y - y- 



a{r) 



where r = Tp -\- y . The normalization constants C and t; are given by 

/2 

and C 



V 



i/ri 



1/u 



(12) 



(13) 



We will see below that the exponential model j/ = 1 provides an adequate fit for both the data and for JV-body models; in 
this case, ri = ^/2 and C = l/^/2. 

The role of ■Ui2('') can be best understood by examining the first momentum moment of the second BBGKY equation, 
which expresses the conservation of particle pairs (Davis & Peebles 1977; Peebles 1980, Equation 71.6), 

1 







(14) 



dt dx 

where a; is a comoving length, a is the scale factor, and Vi2{x,t') is the proper, not comoving, peculiar velocity. 

Let us solve for Vi2 of the dark matter on large scales; the biasing of galaxies enters in a subtle way, as we shall see. In 
the limit ^ < 1, 



Jo 



dx' x'^ ^[x' , t) oc i?^(i) 



(15) 



where i?(i) is the linear growth rate of dark matter fluctuations appropriate for a given cosmology. Equations 14 and 15 yield 



vi2{x,t) = -f{Q)a Ho x^{x) 
where /(H) = din D / din a and 

i{x) = ~S J dx' x'^${x') 



(16) 



(17) 



This equation relates the streaming and correlation function of dark matter, while what is observed are galaxies. The 
quantity Vi2 is a density weighted mean of the pairwise velocities (cf., the discussion in Bertschinger 1992). That is, if expressed 
in terms of a volume-weighted mean: 

^12(7-) = {-[(v(x)-v(x + r))-r] [l + 5(x)][l + 5(x + r)]>,oi„^e . (18) 

We can expand this into a term independent of 6, a term proportional to 6, and a term proportional to 6^ . The latter 
term can be dropped in the limit 5 <C 1. In linear theory, the volume- weighted pairwise velocity field is symmetric, and thus 
the zeroth order term drops out. Thus 



^^12(7-) = {[- (v(x) - v(x + r) ■ r)] [6{x) + 6{x + r)]>^oi^ 



(19) 



Because the velocity field v is the same for dark matter and galaxies, in a biasing model in which 6 galaxies = b6dark matter, 
we have V12 galaxies = bvi2 dark matter- Similarly, ^galaxies = b'^ ^dark matter, and thus Equation 16 yields 



(20) 



galaxies 

2 2 

^12 galaxie s(^X^ — ~P(l Ho X ^galaxie s(^X^ , 

where /3 = /(Q)/fe. A slightly more accurate approximation would be to write: 

2 2 

V12 galaxies{x) = -/3a Ho X ^galaxie s{x) / (1 + ^{x)) , 

and indeed, this is what we shall use in this paper. For convenience, in the remainder of this paper we will drop the subscript 
galaxies on ■Ui2('') and ^(t'), bearing in mind that these quantities refer to the galaxy distribution and not the mass. 



(21) 
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normalizing factor, although we will see below that the signal to noise ratio in this statistic drops below unity on scales above 
~ 10/i"^Mpc. Unfortunately, even on scales as large as 15 h '^Mpc, Q{s) is quite sensitive to the value of the dispersion on 
small scales (cf.. Figure 7) and therefore we are forced into a simultaneous fit for both the dispersion and the streaming. 

On small scales, the clustering is presumed to be stable in the sense that the average separation of galaxy pairs is constant 
in physical coordinates; consequently the first moment of the velocity distribution function (streaming motion) approaches 
zero for small pair separations. Thus in order to isolate the effect of small-scale velocity dispersion (the second moment), we 
look at a statistic with most of its weight on small scales. With perfect data, it would be ideal to take a cut through the 
contour plot at Vp = 0, but in practice, this is very noisy. Thus we sum over the first two bins in Vp to define: 

{^(7-p = l,7r)> = 0.5[^(7-p = 0.5/i"^Mpc,7r) + ^(7-p = 1.5/i"^Mpc,7r)] . (6) 

As with (5(s), we are mainly interested in fitting the shape of this function in order to constrain the velocity dispersion; 
we wish our fitting to be insensitive to the amplitude of ^('•p,7r). Therefore, in analogy with the normalization of (5(s), we 
normalize the statistic in Equation 6, 

Jo (i[rp = l,T))dT 

where Tmax is the maximum value of tt used in the analysis; in the results presented below Tmax =15 /i~''Mpc. 

The solid connected points in Figure 2b plot ^(7r) for the full sample, while the open circles and stars are for the Northern 
and Southern portions, respectively. The three sets of points agree with one another very well, implying that this statistic 
is more robust than ^('•p,7r) as a whole. The quantity ^(7r) is almost a straight line in this linear-log plot, and is positive to 
~ 10/i~''Mpc, beyond which the errors (which we will quantify below) dominate. 



3 MODELING ^(7-p,7r) 

In Paper 1, we projected the contour maps of ^(7'j,,7r) onto the Vp axis to extract the underlying real space correlation function 
^(7-) (Equation 4). However, we can extract more information from our plot of ${rp, tt) than ^(r) alone, as the distortions in 
redshift space contain information on the velocity distribution function of pairs of galaxies. Let f (w|r) be the distribution 
function of velocity differences w for pairs of galaxies separated by vector distance r. A convenient model for ^('•p,7r) is 
(Peebles 1980, Equation 76.6), 

l + arp,'K)= [ d'wF{w\r)[l+^{r)] , (8) 



where = + (tt — Wi/HoY t Ho is the Hubble constant, and is the component of w along the line of sight. Thus 
^(7'j,,7r) may be thought of as a convolution of the real space correlation function with the velocity distribution function. This 
expression does not take into account the correlations one expects between particular features in the density and velocity 
fields (e.g., coherent infall into overdense regions). However, because ^('•p, tt) is computed by averaging over all environments, 
such effects should average out to the extent to which our survey volume represents a "fair" sample of the universe. 

A good approximation to Equation 8 in the regime of slowly varying velocity dispersion is (Peebles 1980, Equation 76.12) 

l + i{rp,i,) = j dw,}{w,\r)[l+i\{rl + {i,-w,IHofy'^\] , (9) 

where /(tosIt") is the velocity distribution function averaged over directions perpendicular to the line of sight, 

f{wi\T) = J dwidw2F{w\T) . (10) 

If we wish to fit Equation 9 to the data, we need a model for the distribution function f^wslr). Let us start by pa- 
rameterizing f (w|r) in terms of its first and second moments; as we will see below, it is these moments which will allow 
us to draw conclusions about Q. The first moment, the mean relative velocity of galaxy pairs of separation r, we label 
■Ui2(r) = — {[v(x) — v(x + r)] ■ r), where the average is number weighted, not volume-weighted. By symmetry this mean 
streaming must be directed along the separation vector of the pair, and is thus a function only of the magnitude of the vector 
r; we shall adopt a convention that assigns positive streaming for infall, in opposition to Hubble expansion. The corresponding 
first moment of / is {^3} = y vi^ij) I t , where y is the component of r in the line-of-sight direction. The second moment of 
F is the 3-d relative velocity dispersion of pairs, £(?•) = {(v(x) — v(x + r) + vi-^ij^i'fY^'^. In general, we should distinguish 
between the longitudinal and transverse velocity dispersion (cf., § 4.2); for now, however, we assume an isotropic velocity 
dispersion tensor. Under this approximation of isotropy, the second moment of / is just the 1-d dispersion, f^(7') = S^(7')/3. 
These moments figure prominently in the kinetic theory formulation of clustering dynamics (Davis & Peebles 1977; Peebles 
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Figure 2. a). The normalized quadrupole moment Q{s^ for tlie full IRAS sample (filled circles), the Northern (open circles) and Southern 
(stars) subsamples, and the full sky with clusters boosted (triangles). Note that Q{s) becomes negative for s > 4/i~^Mpc. b). Same as 
in a) for ^(tt). 

Hamilton (1992) has emphasized the use of multipole moments of the ^('•p,7r) maps to measure the streaming on large 
scales. We shall here limit ourselves to consideration of the quadrupole, which is the major distortion expected as we will see 
below (Equation 23). Consider the statistic (5(s) defined by 

1 

Q{s) ^ 5 ^ , (5) 



where /i = -k/s, and Vni^fi) is the Legendre polynomial of order n. The (5(s) defined here is the quantity ^2(5) defined by 
Hamilton (1992), normalized by the zeroth moment. We apply this normalization in order to minimize sensitivity to the 
amplitude of ^(t-j,, tt), and thus maximize the sensitivity to the redshift distortions themselves. 

Figure 2a shows the (5(s) statistic, for the full sky analysis (solid connected points), the North and South subsamples 
(open circles and stars, respectively), and the boosted sample (triangles). Note that (5(s) is positive for small s, but becomes 
negative, corresponding to an oblate distortion of ^(rj,,7r) along the Tp axis, for s > 4/i~''Mpc. Moreover, the (5(s) for the full 
sky, the two subsamples, and the boosted sky agree quite well. (5(s) does not asymptote to zero at large scales, because of the 
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b) Full Sample-Boosted 
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Figure 1. The IRAS 1.2 Jy ^(rp,7r) for the a) full sample, b) the full sample with clusters boosted, c) the Northern subsample (6 > 5°), 
and d) the Southern subsample (b < 5°). All galaxies in the sample with redshifts between 500 and 30,000 km s"-"- are used. The contours 
are spaced at intervals of 0.1 dex above ^ = 1 (the heavy solid curve); for smaller values of ^, the contours are spaced at arithmetic 
intervals of 0.1. The ^ = contour is the heavy dashed curve, and negative contours are light dashes. The behavior of the redshift space 
correlation function, ^(s), is represented by dotted circular arcs with radii equal to the values of s at which ^(s) = 1.0,0.5, and 0.25. 
For graphical clarity, the map has been twice smoothed by a 1-2-1 boxcar in each direction. Note that on small scales, the contours are 
distorted in the tt direction, while on larger scales, a distortion in the rp direction is evident. 
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then define separations which are parallel (tt) and perpendicular [vp) to the line of sight, 

s ■ 1 ,„x 

- = ^ , (2) 

2 2 
Tp = S ■ S — TT 

With these new variables, we can write down a generalized version of the estimator for the correlation function, 

. NDD{Tp,Tr)NRR{Tp,Tr) 

^^^-'^^ = ' ' 

where iVz)Z)(7'p, tt), iVz)ji(7'p, tt), and Nrr^Tpj-k) refer to data-data, data-random, and random-random pairs with separations 
TT and Tp respectively. This estimator was suggested by Hamilton (19936) as less sensitive to uncertainties in the mean density 
of the sample than is the traditional estimator used by Davis & Peebles (1983) and in Paper 1; in practice we obtained almost 
indistinguishable results using the two estimators. Minimum- variance weighting of the pairs is done, as described in Paper 1, 
with an effective depth between 50 and 110/i~''Mpc, depending on separation. 

Figure la shows ^('•p, tt) for the full 1.2 Jy IRAS sample. The ^('•p,7r) in Figure la was computed using all the galaxies 
in the flux-limited sample in the reds hiftt range 500 < cz < 30, 000 km s ^ , with bin widths and separations of 1 h '"Mpc. 
The contours are spaced at intervals of 0.1 dex above ^ = 1 (the heavy solid curve); for smaller values of ^, the contours are 
spaced at arithmetic intervals of 0.1. The ^ = contour is the heavy dashed curve, and negative contours are light dashes. 
The behavior of the redshift space correlation function, ^(s), is represented by dotted circular arcs with radii equal to the 
values of s at which ^(s) = 1.0, 0.5, and 0.25. For graphical clarity, the map has been twice smoothed by a 1-2-1 boxcar in each 
direction (although the model fitting below is done to the unsmoothed ^('■p,'!')). The small scale distortion discussed above is 
evident as the stretching of the ^{tp,t) contours along the tt direction for values of T-p ^ 2/i~''Mpc. The weak compression of 
the contours along the tt axis for Tp ^ 5/i~''Mpc is the signature of the large scale redshift distortion. 

In Paper 1, we projected ^('•p,7r) onto the Tp axis to estimate the real space correlation function; modeling it with a 
power law yielded 



^('') I „ r,e+o.20 7 _ni/r ) ' 
V3.76Io 23 h iMpcy 

on scales r ^ 20/i~''Mpc. We will use this result extensively below in the modeling of ^(7'p,7r) itself. 

Because the correlation function is a pair-weighted statistic, one expects it to be quite sensitive to the nature of the 
highest-density regions in the sample. Infrared-selected galaxies are known to give lower estimates of the galaxy density in 
cluster cores than do optically selected galaxies (Strauss et al. 1992a); as the cluster cores are regions with the highest pairwise 
velocity dispersion, this may cause us to underestimate the velocity dispersion. In order to address this, we calculated ^('•p,7r) 
for a "boosted" IRAS sample, in which galaxies in the cores of the seven clusters in Table 2 of Strauss et al. are given double 
weight. The results are shown in Figure 16. Reassuringly, this figure differs little in Figure la; we will see that the two yield 
very similar estimates for the velocity dispersion on small scales (§5). 

Sampling fluctuations are a major source of noise in the estimate of ^('•p,7r); there are a finite number of filaments 
and superclusters within the 1.2 Jy survey volume, with a finite number of orientations. Thus note the peaks in ^('•p,7r) 
along the tt axis at tt Ri 12,17, and 22/i~''Mpc, and the extension of the ^ = 0.3,0.4, and 0.5 contours in the tt direction 
at Tp Ri 8/i~''Mpc, which presumably are due to particular structures in the sample which have not been averaged out. To 
gauge the reproducibility of the features of ^(7'p,7r), we subdivide the IRAS sample into a North and South catalogue, split at 
Galactic latitude 6 = 0. Figures Ic and Id show ^('•p,7r) for these two subcatalogues. Roughly speaking, the two subsamples 
are very similar for ^('•p,7r) > 0.4; the discrepancies at smaller amplitudes are indicative of the noise of the derived ^('•p,7r) 
on large scales. The ^('•p,7r) = 1 contour is considerably more prolate in the Northern than Southern catalogue, presumably 
due to the closer proximity of the clusters within the Northern hemisphere. Indeed, our model fits below indicate that the 
velocity dispersion generating this "Finger of God" is larger in the North than the South. 



2.2 Projections of the ^(t"p, tt) Plot 

We will find it useful in fitting models to the data to work with two projections of the ${rp, tt) contours, which are less affected 
by the finite volume of the survey than ^('•p,7r) itself. Our primary aim will be to constrain the first and second moments 
of the velocity distribution function, and thus we choose projections of ^(7'p,7r) which have maximum sensitivity to each of 
these. 

t All redshifts are corrected for the motion of the sun with respect to the barycenter of the Local Group, following Yahil, Tammann, & 
Sandage (1977); no further corrections are applied. 
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1 INTRODUCTION 

Redshift surveys of galaxies can provide statistical estimates of the departures from Hubble flow without the use of direct dis- 
tance measurements of individual galaxies. The prominent elongations in the redshift distribution of cluster galaxies ("Fingers 
of God") are the most conspicuous examples of these "peculiar velocities", but the effect is detectable even in the absence of 
prominent clusters. In particular, the correlation function measured in redshift space ^(s) differs from that in real space ^(r) 
due to two effects: on small scales, ^(s) is suppressed due to the elongation of dense structures along the line of sight in redshift 
space, as mentioned above, while on large scales, coherent infall into overdense regions and outflow out of underdense regions 
enhances the correlations (Kaiser 1987; Fisher et al. 1993a). These effects can be seen in redshift survey data by plotting the 
correlations as a function of the components of the separation of a pair of galaxies in the plane of the sky [vp) and along 
the line of sight (tt). Statistically significant deviations of ^ from concentric circles in the (?•;;, tt) plane are due to redshift 
space distortions: on small scales, we expect the contours to be elongated along the tt axis (a prolate distortion, considering 
^(7'j,,7r) to have rotational symmetry about the tt axis), while on large scales, coherent infall into overdense regions elongates 
the contours along the Vp axis, giving an oblate distortion. 

Previous analyses of redshift space distortions from survey data (Rivolo & Yahil 1981; Davis & Peebles 1983; Bean et 
al. 1983; de Lapparent et al. 1988; Hale-Sutton et al. 1989; Mo et al. 1993) have been primarily concerned with the effect of 
virialized motions of close pairs. On larger scales, a number of authors have made linear theory predictions of the enhancement 
of the correlation function (Sargent & Turner 1977; Kaiser 1987; Lilje & Efstathiou 1989; McGill 1990; Hamilton 1992), which 
have been tested to some degree in JV-body models (Lilje & Efstathiou 1989; Gramann, Cen, & Bahcall 1993). In particular, 
linear theory predicts that the angle-averaged correlation function ^(s) measured in redshift space is enhanced over the real 
space correlation function, ^(t'), by a factor 

if(/3)= [l + ^/3+i/3=] , (1) 

where /3 is related to the linear theory growth factor i?(i) by /3 = din D/dln a Ri Q°'^/fe, and b is the linear bias parameter 
(assumed throughout this paper to be independent of scale) which relates the fluctuations of the galaxy counts to those in 
the underlying matter distribution, [An/n) = b[6p/p). Hamilton (1993a) has detected the large scale infall using angular 
moments of ^('•p,7r) in a redshift survey of galaxies detected by the Infrared Astronomical Satellite [IRAS), flux-limited to 
1.936 Jy at 60/im (Strauss et al. 19926); he concludes that /3 = 0.66^q'22. 

In a previous paper (Fisher et al. 19936; hereafter Paper 1), we presented measurements of the amplitude of the correlation 
function ^('•p,7r) for a sample of 5313 IRAS galaxies flux-limited to 1.2 Jy at 60/im, selected from the IRAS database (Strauss 
et al. 1990; Fisher 1992). The sample covers 87.6% of the sky, and probes to an effective depth of Ri 8000 km . In this paper, 
we focus on the redshift space distortions in ^('•p,7r), in order to characterize the amplitude of the velocity field of galaxies. 
The quantity ${rp, tt) is related to ^(r) through the distribution function of the difference between the peculiar velocity of 
pairs of galaxies at a given separation, f (w|r). This quantity is of fundamental interest, because on large scales, linear theory 
relates the first moment of F to the cosmological density parameter Q and the bias parameter b. On small scales, the Cosmic 
Virial Theorem connects the second moment of f to Q and b. Thus in § 2, after presenting the IRAS ^('•p,7r) itself, we 
define two projections of ^(7'j,,7r) designed to be maximally sensitive to the two moments; we will fit models directly to these 
projections. 

In § 3, we present the formalism needed for modeling ^('•p,7r), and show the predictions of linear theory and the Cosmic 
Virial Theorem. Both the shape of F and the scaling of its moments with separation are unknown. Rather than trying to 
extract these unknown functions from the data, we turn to JV-body simulations of several variants of Cold Dark Matter 
(CDM) models to characterize their forms (§ 4.1). In addition, we extract many Monte-Carlo realizations of the observational 
sample from one of the JV-body simulations in order to characterize the errors in the observed quantities, thus allowing a 
fit of models for the redshift distortions to the data (§ 4.2). We minimize with respect to the unknown form of f at a 
given separation, as well as the normalizations of the first and second moments of the velocity distribution function. We test 
our methods with the Monte-Carlo simulations, and show that we recover unbiased estimates of the normalizations of the 
velocity moments (§ 4.3). In § 5, we apply these techniques to the IRAS data themselves, and discuss their implications for 
Q and in § 6. We conclude in 7. 



2 THE OBSERVED QUANTITIES 
2.1 The Observed ^(T-p,7r) 

As described in Paper 1, we quantify the effects of the redshift space distortions by computing the correlation function as a 
function of separations parallel and perpendicular to the observer's line of sight. Explicitly, given a pair of galaxies with redshift 
positions, Vi and V2, we define a separation in redshift space, s = Vi — V2, and an observer's line of sight, 1 = ^{vi +V2). We 
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ABSTRACT 

We examine the effect of redshift space distortions on the galaxy two-point correlation 
function ^(r"p,7r) as a function of separations parallel (r-p) and perpendicular (tt) to 
the line of sight. Modeling ^(r"p,7r) measured from a full-sky redshift survey of IRAS 
galaxies allows us to characterize the moments of the velocity distribution function 
of pairs of galaxies. We are guided in our parameterization of models by results from 
numerical simulations of Cold Dark Matter (CDM) models. 

It is essential that the model for (.{vp, tt) contain the effects of both the first and second 
moments of the velocity distribution function, as they distort the redshift space cor- 
relations in opposing directions. We develop a method of fitting for the amplitudes of 
the velocity moments, and show that we can recover the correct values in Monte-Carlo 
realizations of the data drawn from TV-body simulations. We find that the relative ve- 
locity dispersion of pairs of IRAS galaxies is cr(T") = SlT^^g km at t" = 1 h~^M.^c, 
consistent with previous estimates derived from optically selected galaxy catalogues. 
Unfortunately, the use of this result to estimate fi via the Cosmic Virial Theorem is 
thwarted by large systematic uncertainties, making the application of this theorem to 
existing redshift surveys of little value. 

We also fit for the mean relative streaming velocity of pairs, ■Ui2(f"), which describes the 
growth of fluctuations on both linear and nonlinear scales. We find that ■Ui2(f") = IBT^gy 
km at T" = 4/i~"'^Mpc, so that on average, approximately half the Hubble expansion 
velocity of pairs at this separation is canceled by infall. At t" = 10 h~^M.^c, the amplitude 
of the streaming is lower and ■Ui2(f") = km s"-*^. Linear perturbation theory then 

implies that Q°-^/h = QAhto^^s on scales ~ 10 - IS/i'^Mpc. The amplitude of Vi2{r) 
is sensitive to the assumed shape of <j{r); if the latter deviates substantially from a 
virialized form on small scale, our best fit amplitude of Vi2{f) can deviate by a factor of 
two. Our derived result for /3 is intermediate between that found on ~1 h~^M.^c scales 
and on 15 — 40 /i~"'^Mpc scales, arguing either for a strong dependence of /3 on scale, or 
an error in the determination of (3 on one or more scales. 

Key words: Cosmology: large-scale structure 
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